Exponential peak and scaling of work fluctuations in modulated systems 
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We extend the stationary-state work fluctuation theorem to periodically modulated nonlinear 
systems. Such systems often have coexisting stable periodic states. We show that work fluctuations 
sharply increase near a kinetic phase transition where the state populations are close to each other. 
The work variance is proportional here to the reciprocal rate of interstate switching. We also show 
that the variance displays scaling with the distance to a bifurcation point and find the critical 
exponent for a saddle-node bifurcation. 
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Since their discovery in the early 90s [H, 0, Q, fluc- 
tuation theorems have been attracting increasing inter- 
est. They establish general features of fluctuatingsys- 
tems away from thermal equilibrium, see Refs. 0,13 for 
reviews. A major "test bed" for fluctuation theorems 
is provided by dynamical systems with a few degrees of 
freedom coupled to a thermal bath, a Brownian particle 
being an example. Much of the corresponding theoreti- 
cal and experimental work refers to (i) modulated linear 
systems, where fluctuations have been studied both in 
transient and stationary regimes [1, 0, H, S [13, EIlj , 
and (ii) nonlinear systems, initially at thermal equilib- 
rium, driven to a different, generally nonequilibrium state 

Fluctuations in nonequilibrium dynamical systems 
have been attracting attention also in a different con- 
text. They play an important role in various types of 
mesoscopic vibrational systems of current interest. Be- 
cause damping of the vibrations is typically weak, even 
a moderately strong resonant force can excite them to 
comparatively large amplitudes, where the nonlinearity 
becomes substantial. As a result, the system may have 
two or more coexisting stable states of forced vibrations 
[l^ . Fluctuations can cause switching between these 
states [l^ and thus significantly affect the overall be- 
havior of the system even where they are small on aver- 
age. Many features of the switching behavior and a range 
of phenomena and applications related to the switching, 
from quantum measurements to resonant frequency mix- 
ing and to high-frequency stochastic resonance have been 
studied experimentally [13, [U, [H, [H H [H, [H, [H [H . 

In this paper we analyze work fluctuations in periodi- 
cally modulated nonlinear dynamical systems coupled to 
a bath. We derive the stationary state work fluctuation 
theorem and show that, under fairly general assumptions, 
the distribution of fluctuations of work done by the mod- 
ulating force over a long time t is Gaussian. In common 
with systems close to thermal equilibrium, the work vari- 
ance (T^ is proportional to the average work (W), but the 
proportionality coefficient is not universal and depends 
on system parameters. It becomes exponentially large in 
bistable systems in the range of a kinetic phase transi- 
tion where the stationary populations of the vibrational 
states are close to each other. This parameter range has 



similarity with the region of a first-order phase transition 
where molar fractions of the coexisting phases are close 
to each other pj3, 29, 3l|. 

Large work fluctuations are related to the difference 
in the power absorbed from the force in different vibra- 
tional states i = 1,2. Switching back and forth be- 
tween the states leads to significant fluctuations of the 
absorbed power. Their correlation time is determined by 
the switching (transition) rates. For a small character- 
istic intensity D of the noise that comes from the bath, 
these rates are small compared to the dynamical relax- 
ation rate in the absence of noise and the modula- 
tion frequency ajp. Then of primary interest are period- 
averaged transition rates Vij. They often display activa- 
tion dependence on D, with i>ij oc exp{—Ri/D), where 
Ri is the characteristic activation energy of a transition 
i j- Since work fluctuations accumulate power fluctu- 
ations, and the typical accumulation time for interstate 
fluctuations is ~ ^ij^ j the exponential smallness of Vij 
leads to an exponentially large factor in the variance a^. 

We present the results for a most simple model: a non- 
linear classical dynamical system modulated by a peri- 
odic force F{t) = ^„ F(n) exp(mw_Ft); the coupfing en- 
ergy is ~F{t)q, where q is the system coordinate. The 
system is additionally coupled to a bath, which leads to 
relaxation and fluctuations. Work done by the force over 
time T is 



W = W{t) = / dtF{t)q{t) 



(1) 



We are interested primarily in steady-state fluctuations, 
i.e., we assume that the system had come to the steady 
state well before the time t = when the work ([1]) started 
to be measured. This steady state is periodic in time 
with modulation period Tp = 21^ /top. We further assume 
that correlations in the system decay sufficiently fast, for 
example, exponentially for long times, and that the time 
T largely exceeds the characteristic correlation time of 
the system tcorr- Often for bistable systems ^corr ^ ^/^ij- 
Work fluctuations can be expressed in terms of the 
correlation function of velocity fluctuations Q{t,t') = 
{dq{t)dq(t')), where (...) means ensemble average and 
^Qit) ~ <i{t) — (9(0} • Because the system is in a steady 
periodic state, we have Q{t, t') = Q{t + Tp, t' + Tp), and 
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therefore 

Q{t, t') = ^ Q{n; t ~ t') exp [inajpit + t')/2] . (2) 

n 

We first consider the variance of the work distribution 
cr^ = cr2(T-) = {{SW)^), where VK(r) - {W{t)). In 

the Hmit of large r 



n + m 
m — n; — - — ivp 



Q{n;Lu) = (27r)-i / e''"*Q„(i). 



(3) 



Here we have taken into account that the correlation 
functions Q{n]t — t') decay on time icon- much smaller 
than T. Therefore the limits of integration over t ~ t' 
could be extended from — oo to oo. 

Decay of correlations on a time scale small com- 
pared to T allows one to simplify the expressions for 
higher-order moments of 6W in a standard way. The 
3rd moment {{SW)^) is determined by an integral over 
ti,t2,t3 of the correlator {Sq{ti)Sq{t2)Sq{t3)) . Because 
^1 — ^21,1*1 — < ^corr, we havc {{SW)^) oc T for 
large r, and therefore {{5Wf)/{{SW)Y^^ « t'^/^, 
i.e., the 3rd moment is small for large r. The 4th 
moment {{6W)^) is determined by the integral of the 
correlator {6q{ti)Sq{t2)Sq{t3)Sq{t4)) . The main contri- 
bution to this integral comes from decoupling the cor- 
relator into pairs {6q{tnj^)Sq{tn.2)){Sq{t„^)5q{tnJ) with 



^n2 1 1 



tr 



while It 



in3 I 



1, 



,4). This gives {{SWf) 



rii 



n.3 I 

2)2 OC t2. 

Higher-order correlations in {{5W)'^) give a compara- 
tively small contribution oc t. The analysis can be im- 
mediately extended to higher moments of 5W . It shows 
that the overall distribution of work fluctuations P{W) 
is Gaussian, 

P{W) ^ (27r(T2)-i/2exp[-(W^-(W^))2/2a2] , (4) 

It follows from Eq. © that P{W)/P{-W) = 
e'XY){2W {W) / a^) ^ as in the stationary state fluctuation 
theorem for systems close to thermal equilibrium and for 
modulated linear systems, and the variance of the work 
distribution a"^ (x t o< 



Here, qj{n) is the Fourier component of the coordinate 
qj{t) in a periodic state j, qj{t) = iji'n) e'xp{inijj pt) . 

In contrast to the average work, the variance a'^ has 
contributions of two different types. One comes from 
small-amplitude fluctuations about the stable states. It 
is given by the sum of partial variances a\ 2 weighted 
with the state populations. The variances a\ 2 can be 
obtained by linearizing equations of motion about the 
corresponding stable state and can be written as 



a2 = QW,, J = 1,2. 



(6) 



The constant Cj depends on the details of the system 
dynamics. For a system coupled to a thermal bath at 
temperature T, in the weak-field limit we have a'^ = 
2kBTWj. However, for strong field this relation does 
not hold in nonlinear systems, generally. 

Fluctuations about qj become large near a saddle-node 
bifurcation point where the state j disappears. Here, one 
of the motions of the system is slow [3l|, i.e., there is a 
"soft mode" . Near the bifurcational (critical) position of 

the periodic state qj''\t) one can quite generally write 

q{t) - qf'\t) = qsm{t)>c{t), where M:{t) = M:{t + tf) is a 
periodic function, whereas qsui{t) is a slowly varying am- 
plitude that depends on initial conditions. After rescaling 
the equation of motion for qsm{t) can be written as [32| 



V + f{t), {f{t)f)t') = 2D'5{t~t'). (7) 



Here, rj is the distance to the bifurcation point, for exam- 
ple, the scaled difference between the amplitude A of the 
field and its critical value A,, at the bifurcation point. 
The noise f{t) is assumed to be white because of the 
slowness of gsm(t); its intensity is D' oc D. 

For ?7 > in the absence of noise the system ([7]) has 



a stable state where Qs 



^/2. Small fluctuations 



about this state have variance D'rj ^/2/2 and decay over 
time = 7/-1/2/2. Therefore, from Eq. ©, near the 
bifurcation point where a state j disappears 



(8) 



{W). However, for strong pe- p^^^^^ independent of D and 77, and W, does not 



riodic modulation there is no known general expression 
that would relate the average velocity {q{t)) to the mod- 
ulating force in terms of the correlation functions Qn{'-o)- 
Moreover, as we show, the ratio a"^ / (W) may display 
sharp narrow peaks as a function of system parameters. 

For weak noise, a system with two coexisting sta- 
ble states j = 1,2 mostly performs small fluctuations 
about these states and only occasionally switches be- 
tween them. Then to leading order in the noise intensity 
D the average work is a sum of partial works Wi_2 in each 
of the states weighted with the stationary populations of 
the states z/;|*2 [19| . 

{w) = ' ^1 = E 

J = l,2 n 

=V2il {v]_2^V2i), w^' = l-u;f. (5) 



diverge for r/ ^ 0. 

Equation ([5]) shows that the partial work variance 
scales as if- with the distance to the bifurcation point. 
The critical exponent is ^ = — 1. 

The other contribution to comes from fluctuation- 
induced transitions between the states. The transitions 
lead to fluctuations of the state populations Wj (t) . These 
fluctuations are slow, 



{6wi{t)dwi{t')) 



wfw2 exp [-Z^trl* - t'\] 
^12 + V21, 



(9) 



where 5wi{t) = wi{t) — wf = —Sw2{t). In turn, they lead 
to slow fluctuations of the velocity q{t) w ^ . qj{t)wj{t) 
with decay time given by the reciprocal total transition 
rate f^"^ 3> tf- 
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From Eqs. (P, the contribution to the 

work variance from interstate transitions is 

al « Mii^trT)-^ (Wi - W2f , M = 2«, (10) 
and the total variance is 

-' = E, + (11) 

Equation (fTU)l is the central result of the paper. It 
shows that the work variance is proportional to the 
squared difference of the partial works in the stable states 
and is inversely proportional to the transition rate i^tr- 
The rate ftr cx exp (— min^ Ri/D) is exponentially small 
for small noise intensity. Respectively, the variance pU]) 
can be exponentially large compared to the variance due 
to small fluctuations about attractors ©. 

Factor A4 in Eq. pO)) sharply depends on the param- 
eters of the system and the field F{t). It is small, 

M (xexp[-\Ri- R2\/D], (12) 

except for a narrow range of the kinetic phase transi- 
tion where the transition activation energies are close to 
each other, \Ri — R2\ ^ D. At its maximum Ad = 1/2. 
Factor A4 determines also the intensity of super nar- 
row peaks (of width i^tr ^7^) the power spectra 
of modulated bistable systems and the spectra of ab- 
sorption/amplification of an additional field [l^, [s^l- Its 
exponential dependence on the distance to the kinetic 
phase transition was seen in the simulations [s^ . [ssj and 
experiment [2^, [s^ . 

The average work (W), Eq. switches between the 
values Wi and W2 when the system goes through the 
kinetic phase transition. Therefore the ratio a'^/{W) cx 
M.VtT:^ displays a sharp peak at the transition. 

The above results can be illustrated using as an exam- 
ple a resonantly driven underdamped Duffing oscillator, 
a model that describes a number of systems studied in 
recent experiments. In the absence of noise the oscillator 
dynamics is described by the equation 

q + ulq + ^q^ + 2Tq = AconLOFt. (13) 

We assume that the detuning of the field frequency from 
the oscillator eigenfrequency 5uj = ljf — ojq and the 
friction coefficient F are small: |(5a;|,F <C luq. Then 
the oscillator can display bistability of forced vibra- 
tions already for a comparatively small driving ampli- 
tude A, where the vibrations remain almost sinusoidal, 
qj{t) K. ttj cos{uJFt + {j ^ 1,2). ExpHcit expressions 
for the amplitudes ai^2 and phases (j}i^2 are well known 
[r§ |. The interstate transition rates are well understood 
also [ll,[2l|,|23,[2l,p. 

The partial work in a stable vibrational state j is Wj ~ 
Frwl^a^. By calculating the partial variances a| due to 
thermal noise one obtains 

-g- = 2kBTZr^ [{Z, - 2f + 4^(2 (y, - 1)^] , (14) 

n = 5uj/T, = l + fi2(Y,- - l)(3Yj- - 1), 



where Yj = 3'ya'j /SlofSuj is the scaled squared vibration 
amplitude. Eq. (fH)) refers to the case where fluctua- 
tions and friction come from coupling to the same ther- 
mal reservoir; in fact, it is not limited to the model (jl3p 
and applies in a general case where the density of states 
of the reservoir weighted with the interaction is smooth 
near ujf, cf. [Tgj . 

It follows from Eq. (fTi|) that in the linear-response 
regime, where Yj oc A"^, Yj ^ 1, the ratio cr'j/Wj 
2kBT is given by the standard stationary state work 
fluctuation theorem. This is in agreement with the re- 
cent theoretical and experimental results on a periodi- 
cally modulated linear system (llj . However, for stronger 
driving the ratio is no longer given by 2A:bT. When 
the driving amplitude A or frequency lof approach their 
bifurcational value, the ratio diverges, 

a^/W,^2kBTGj{n)ir\ v^P-f3^i\n). (15) 

Here, (3 = SjA"^ /32ojp{6Lu)^ is the reduced field ampli- 
tude, Pb '^^ are the bifurcational values of /3 [l8[, and 

Gj .= n-^p'^^^/YjBi^YjB - 2) (Y,B is the bifurcational 
value of the reduced squared vibration amplitude). In 
agreement with Eq. ([5]), (Jj /Wj scales as r]~^ with the 
distance 77 to the bifurcation point. The full dependence 
of a'^/Wj on (3 is shown in Fig. [TJ 
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FIG. 1: (Color online) Scaled ratios of the partial work vari- 
ance to mean partial work Sj = (Tj/2fcsTWj, j — 1,2, as 
functions of the reduced squared modulation amplitude /3 for 
a Duffing oscillator; = F/Slo — 0.3. The curves 1 and 

2 refer to large- and small-amplitude vibrations, respectively. 
The functions 51,2 diverge at the corresponding bifurcation 
points of the oscillator, /3^^' ^ 0.088 and p'g'' a; 0.18. 

It should be possible to see the scaling (fT5|) for 
modulated oscillators for the same conditions as the 
scaling of the switching activation energies P 111. 
Similarly, the exponential peak of the work variance 
should be seen in experiments analogous to those where 
there were observed other kinetic phase transition phe- 
nomena like super narrow peaks in the power spec- 
tra, high-frequency stochastic resonance, and fluctuation- 
enhanced frequency mixing [25l . [28l . [36| . 
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Experiments on the scaling and on the kinetic phase 
transition should be conducted in a different way. The 
scaling can be observed in the quasistationary regime. 
The system should be prepared in the metastable state 
near a bifurcation point, and the duration of the exper- 
iment should be shorter than the lifetime of the state. 
In contrast, the exponential peak of the variance can be 
seen provided the duration of the measurement exceeds 
the reciprocal switching rate v^^^. 

In conclusion, we have considered the stationary-state 
work fluctuation theorem for dynamical systems modu- 
lated by a strong periodic field. We have shown that, 
if the system has coexisting states of forced vibrations. 



work fluctuations may become strong. The ratio of their 
variance to average work is proportional to the reciprocal 
rate of interstate transitions. It displays a sharp expo- 
nentially high peak as a function of the distance to the 
kinetic phase transition. In a different parameter range, 
near a saddle-node bifurcation point where one of the sta- 
ble vibrational states disappears, in the quasistationary 
regime the variance displays scaling with the distance to 
the bifurcation point. 
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stimulating discussion of the fluctuation theorems. This 
work was supported in part by the NSF through grant 
PHY-0555346. 



[1] D. J. Evans, E. G. D. Cohen, and G. P. Morriss, Phys. 

Rev. Lett. 71, 2401 (1993). 
[2] D. J. Evans and D. J. Searles, Phys. Rev. E 50, 1645 

(1994). 

[3] G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 
2694 (1995). 

[4] D. J. Evans and D. J. Searles, Adv. Phys. 51, 1529 
(2002). 

[5] C. Bustamante, J. Liphardt, and F. Ritort, Phys. Today 

58 (7), 43 (2005). 
[6] G. M. Wang, E. M. Sevick, E. Mittag, D. J. Searles, and 

D. J. Evans, Phys. Rev. Lett. 89, 050601 (2002). 
[7] R. van Zon, S. Ciliberto, and E. G. D. Cohen, Phys. Rev. 

Lett. 92, 130601 (2004). 
[8] N. Garnier and S. Ciliberto, Phys. Rev. E 71, 060101 

(2005). 

[9] F. Zamponi, F. Bonetto, L. F. Cugliandolo, and J. Kur- 
chan, J. Stat. Mech. p. P09013 (2005). 
[10] E. G. D. Cohen and R. van Zon, cond-mat/0701683 
(2007). 

[11] S. Joubaud, N. Garnier, and S. Ciliberto, cond- 

mat/0703798 (2007). 
[12] A. Imparato, L. Peliti, G. Pesce, G. Rusciano, and 

A. Sasso, cond-mat/0707.0439 (2007). 
[13] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997). 
[14] G. E. Crooks, Phys. Rev. E 60, 2721 (1999). 
[15] G. Hummer and A. Szabo, Proc. Nat. Acad. 98, 3658 

(2001). 

[16] D. Collin, F. Ritort, C. Jarzynski, S. B. Smith, I. Tinoco, 

and C. Bustamante, Nature 437, 231 (2005). 
[17] V. Blickle, T. Speck, L. Helden, U. Seifert, and 

C. Bechinger, Phys. Rev. Lett. 96, 070603 (2006). 
[18] L. D. Landau and E. M. Lifshitz, Mechanics (Elsevier, 

Amsterdam, 2004), 3rd ed. 
[19] M. I. Dykman and M. A. Krivoglaz, Zh. Eksp. Teor. Fiz. 

77, 60 (1979). 

[20] L. J. Lapidus, D. Enzer, and G. Gabrielse, Phys. Rev. 
Lett. 83, 899 (1999). 



[21] I. Siddiqi, R. Vijay, F. Pierre, C. M. Wilson, L. Frun- 
zio, M. Metcalfe, C. Rigetti, and M. H. Devoret, cond- 
mat/0507248 (2005). 

[22] J. S. Aldridge and A. N. Cleland, Phys. Rev. Lett. 94, 
156403 (2005). 

[23] K. Kim, M. S. Heo, K. H. Lee, H. J. Ha, K. Jang, H. R. 

Noh, and W. Jhe, Phys. Rev. A 72, 053402 (2005). 
[24] C. Stambaugh and H. B. Chan, Phys. Rev. B 73, 172302 

(2006). 

[25] C. Stambaugh and H. B. Chan, Phys. Rev. Lett. 97, 
110602 (2006). 

[26] I. Siddiqi, R. Vijay, M. Metcalfe, E. Boaknin, L. Frunzio, 
R. J. Schoelkopf, and M. H. Devoret, Phys. Rev. B 73, 
054510 (2006). 

[27] B. Abdo, E. Segev, O. Shtempluck, and E. Buks, J. Appl. 

Phys. 101, 083909 (2007). 
[28] R. Almog, S. Zaitsev, O. Shtempluck, and E. Buks, Appl. 

Phys. Lett. 90, 013508 (2007). 
[29] R. Bonifacio and L. A. Lugiato, Phys. Rev. Lett. 40, 1023 

(1978). 

[30] L. A. Lugiato, Prog. Opt. 21, 69 (1984). 

[31] J. Guckenheimer and P. Holmes, Nonlinear Oscillators, 
Dynamical Systems and Bifurcations of Vector Fields 
(Springer- Verlag, New York, 1987). 

[32] M. 1. Dykman, B. Golding, and D. Ryvkine, Phys. Rev. 
Lett. 92, 080602 (2004). 

[33] M. I. Dykman, M. A. Krivoglaz, and S. M. Soskin, 
in Noise in Nonlinear Dynamical Systems, edited by 
P. V. E. McClintock and F. Moss (Cambridge Univer- 
sity Press, Cambridge, 1989), vol. 2, pp. 347-380. 

[34] M. I. Dykman, R. Mannella, P. V. E. McChntock, and 
N. G. Stocks, Phys. Rev. Lett. 65, 48 (1990). 

[35] M. I. Dykman, D. G. Luchinsky, R. Mannella, P. V. E. 
McClintock, N. D. Stein, and N. G. Stocks, Phys. Rev. 
E 49, 1198 (1994). 

[36] H. B. Chan and C. Stambaugh, Phys. Rev. B 73, 224301 
(2006). 



